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Fermionic allialine-earth atoms have unique properties tliat make tliem attractive candidates for tlie 
realization of novel atomic clocks and degenerate quantum gases. At the same time, they are attracting 
considerable theoretical attention in the context of quantum information processing. Here we demon- 
strate that when such atoms are loaded in optical lattices, they can be used as quantum simulators of 
unique many-body phenomena. In particular, we show that the decoupling of the nuclear spin from the 
electronic angular momentum can be used to implement many-body systems with an unprecedented de- 
gree of symmetry, characterized by the SU(N) group with N as large as 10. Moreover, the interplay of 
the nuclear spin with the electronic degree of freedom provided by a stable optically excited state allows 
for the study of spin-orbital physics. Such systems may provide valuable insights into strongly correlated 
physics of transition metal oxides, heavy fermion materials, and spin liquid phases. 



The interest in fermionic alkaline-earth atoms |1|, 

stems from their two key features: (1) the presence of 
a metastable excited state "^Pq coupled to the ground ^5*0 state 
via an ultranarrow doubly-forbidden transition iQl] and (2) the 
almost perfect decoupling lit] of the nuclear spin / from the 
electronic angular momentum J in these two states, since they 
both have J = 0. This decoupling implies that s-wave scat- 
tering lengths involving states ^Sq and ^Pq are independent 
of the nuclear spin, aside from the restrictions imposed by 
fermionic antisymmetry. We show that the resulting SU(N) 
spin symmetry (where iV = 2/ + 1 can be as large as 10) to- 
gether with the possibility of combining (nuclear) spin physics 
with (electronic) orbital physics open up a wide field of ex- 
tremely rich many-body systems with alkaline-earth atoms. 

In what follows, we derive the two-orbital SU(N)- 
symmetric Hubbard model describing alkaline-earth atoms 
in ^5*0 and ^Pq states trapped in an optical lattice. We fo- 
cus on specific parameter regimes characterized by full or 
partial atom localization due to strong atomic interactions, 
where simpler effective spin Hamiltonians can be derived. 
The interplay between orbital and spin degrees of freedom 
in such effective models is a central topic in quantum mag- 
netism and has attracted tremendous interest in the condensed 
matter community. Alkaline earth atoms thus provide, on 
the one hand, a unique opportunity for the implementation 
of some of these models for the first time in a defect-free 
and fully controllable environment. On the other hand, they 
open a new arena to study a wide range of models, many of 
which have not been discussed previously, even theoretically. 
We demonstrate, in particular, how to implement the Kugel- 
Khomskii model studied in the context of transition metal ox- 
ide sjl 1110 [ilB, the Kondo lattice model Qlillil 
lilliltllliallllilllllllllllli] studied in context of 



manganese oxide perovskites 112011 and heavy fermion materi- 
als |25J, as well as various SU(N)-symmetric spin Hamilto- 
nians that are believed to have spin liquid and valence-bond- 
solid ground states iUli IH S lH S S H . For ex- 
ample, we discuss how, by appropriately choosing the initial 
state, a single alkaline-earth atom species with 7 = 9/2 (such 
as ^^Sr) can be used to study experimentally such a distinc- 
tively theoretical object as the phase diagram as a function of 
N for all N < 10. 

Before proceeding, we note that, while an orthogonal sym- 
metry group SO(5) can be realized in alkali atoms IBSll . 
proposals to obtain SU(N>2)-symmetric models with alkali 
atoms ll36l fSTll and solid state systems HH fsstl are a sub- 
stantial idealization due to strong hyperfine coupling and a 
complex solid state environment, respectively. In this context, 
alkaline-earth-like atoms make a truly exceptional system to 
study models with SU(N>2) symmetry. 
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FIG. 1: Interaction parameters between g (green) and e (yel- 
low) atoms loaded in tlie lowest vibrational state of the corre- 
sponding optical lattice. Here we assumed / = 1/2, and the ar- 
rows indicate tlie mi — ±1/2 spin states. \s,t) denote the singlet 
and triplet nuclear spin states of the two atoms (only one of three 
triplet states - | ||) - is shown). The dashed circle represents anti- 
symmetrization of the nuclear spin state (i.e. |s)). The interaction en- 
ergy Ux (X — gg, ee, eg^, eg~) is proportional to the correspond- 
ing scattering length ax ■ 



Information). 



The independence of each of the four scattering lengths 
from the nuclear spin state is essential to the fulfillment of 
the SU(N) symmetry of our model (see next Section). This 
independence is a consequence of the decoupling between 
nuclear and electronic degrees of freedom exhibited during 
the course of a collision involving any combination of g or 
e states, which both have J = 0. While for the |e) = "^Pq 
atom, the decoupling is slightly broken by the admixture with 
higher-lying P states with J 7^ 0, this admixture is very small 
yj] and the resulting nuclear-spin-dependent variation of the 
scattering lengths is also expected to be very small, on the or- 
der of 10^'^ (see Supplementary Information). For Ugg, which 
does not involve state |e), this variation should be even smaller 
(~ 10-9). 



Many-body dynamics of alkaline-earth atoms in an optical 
lattice 



We begin with the Hamiltonian describing cold fermionic 
alkaline-earth atoms in an external trapping potential: 

H = J2l Sv^i^{v){-^V'- + y„(r))*_(r) (1) 

am 

+nojo J dh{p,{r) - p,(r)) + ^i±^ J d^vp,{v)pg{v) 

+ X! Saa j (f r Pamir) p am' (j) 

I d^r^>l^{r)^flAr)^gm'{r)^emir). 
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Here ^'^^(r) is a fermion field operator for atoms in internal 
state jam), where a = g CSq) or e (^Pq) denotes the elec- 
tronic state and m ~ —I, . . . , I denotes one of the N ~ 21+ 1 
nuclear Zeeman states. The density operators are defined as 
/3am(r) = 1'J„„(r)^'a,„(r) andpQ(r) = Y.mParn{r). The 
term V^(r) describes the external trapping potential, which 
we will assume to be an optical lattice independent of the nu- 
clear spin: even for a relatively deep lattice with a 100 kHz 
trap frequency, tensor and vector light shifts should be well 
below 1 Hz III]. hiJct is the transition energy between \g) and 
|e). Extra lasers can be used to drive transitions between \g) 
and |e) levels |[ll|2l. Since we will only need these extra lasers 
for system preparation, we have not included the correspond- 
ing terms in the Hamiltonian. 

The interaction is characterized by four s-wave scattering 
lengths ax, X — gg,ee, eg~^, eg~, which define four interac- 
tion parameters gx = '^nh^ax /M, where M is atomic mass. 



^gg 1 ^ee, 



and afg are the scattering length for two atoms in 



the electronic state \gg), |ee), and |±) = {\ge) + \eg))/V2, 
respectively. As shown in Fig.[T] the fermionic antisymmetry 
then forces the nuclear state to be symmetric for the only anti- 
symmetric electronic state |— ) and antisymmetric otherwise. 
Very few ax are known at the moment (see Supplementary 



The interaction terms in Eq. ([T]i describe the most general 
s-wave two-body interaction consistent with elastic collisions 
as far as the electronic state is concerned and with the inde- 
pendence of the scattering length from the nuclear spin. While 
the assumption of elasticity for g-g and e-g collisions is well 
justified, since no inelastic exit channels exist, e-e collisions 
are likely to be accompanied by large losses, which means that 
the magnitudes of the imaginary and real parts of the e-e scat- 
tering length are likely to be comparable (see Supplementary 
Information). Therefore, we focus below on those situations 
where two e atoms never occupy the same site. 



We assume that only the lowest band in both e and g lat- 
tices is occupied and expand the field operators in terms of 
the corresponding (real) Wannier basis functions ^'Q,„(r) ~ 
J2j WQ,(r — Vj)cjam, where c^j^m creates an atom in internal 
state \am) at site j (centered at position Yj). Eq. ([T]i reduces 
then to a two-orbital single-band Hubbard Hamiltonian 

H = Ja(cL„Cj„,„+h.c.) + ^^nj„(njc«-l) 

{j.i)oL,m j.a 
3 j,rn.m' 

Here Ja = - J d^rwa{r){-j^V'^ + 14(r))u;Q(r - Tq) are 
the tunneling energies, Tq connects two nearest neighbors 
h.c. stands for Hermitian conjugate, rijam 



ja — X] 



m "-jam 



'^lam^jarn, and 

The tunneling is isotropic, which is a cru- 
cial difference between this model and its analogues in solid 
state systems with orbital degeneracy lUb. The sum is 
over pairs of nearest neighbor sites i, j. V = {U^g + U~g)/2 
and Vex ~ [U^g — C^eg)/2 describe the direct and exchange 
interaction terms. The onsite interaction energies are Uaa 



gaa J d^rwair) and Uj; 
terms, proportional to 



rieg 



J d^rwl{r)wg{r). Constant 
, are omitted in Eq. (|2]i. Exper- 
imental control over the parameters in Eq. (|2]i will allow us to 
manipulate the atoms (see Methods). 
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Symmetries of the Hamiltonian 

To understand the properties of the Hamiltonian in Eq. Q, 
we consider its symmetries. We define SU(2) pseudo-spin al- 
gebra via 



j jmaf3 



t M 



(3) 



where cr'* (/i = x, y, z) are Pauli matrices in the {e, 5} basis. 
We further define nuclear-spin permutation operators 



(4) 



r cm 



which satisfy the SU(N) algebra [S"™, S^] 
and thus generate SU(N) rotations of nuclear spins {N = 21+ 
1). 

In addition to the obvious conservation of the total number 
of atoms n = '^jinje+njg), H exhibits U (1) x SU{N) sym- 
metry (see Methods for the discussion of enhanced symme- 
tries), where ?7(1) is associated with the elasticity of collisions 
as far as the electronic state is concerned {[T^, H] = 0) and 
SU(N) is associated with the independence of scattering and 
of the trapping potential from the nuclear spin {[S™', H] = 
for all n, m). The two-orbital SU(N)-symmetric Hubbard 
Hamiltonian in Eq. (|2]i is a generalization to > 2 of its 
SU(2)-symmetric counter part and to two orbitals of its 
single-orbital counterpart 12811 . The SU(N) symmetry and the 
largely independent spin and orbital degrees of freedom are 
two unique features present in alkaline-earths but absent in 
alkalis due to strong hyperfine interactions. 

One important consequence of SU(N) symmetry is the con- 
servation, for any m, of S*™, the total number of atoms with 
nuclear spin m. This means that an atom with large /, e.g. 
^^Sr (/ = 9/2), can reproduce the dynamics of atoms with 
lower I if one takes an initial state with = for some 
m values. To verify SU(N) symmetry of the interaction ex- 
perimentally, one could, thus, put two atoms in one well in 
spins m and m' and confirm that collisions do not populate 
other spin levels. This feature of SU(N) symmetry is in stark 
contrast to the case of weaker SU(2) symmetry, where the de- 
pendence of scattering lengths on the total spin of the two col- 
liding particles allows for scattering into spin states other than 
m and m'. We note that although collisions are governed by 
electronic interactions and obey the nuclear-spin SU(N) sym- 
metry, the nuclear spins still indirectly control the collisions 
via fermionic statistics and give rise to effective spin-orbital 
and spin-spin interactions. 

One can alternatively implement the two-orbital Hubbard 
model with two ground-state species of alkaline-earth atoms 
(e.g. Yb and i^^yb, or ^^^yb and ^^Sr). If we still re- 
fer to them as \g) and |e), the nuclear distinguishability and 
the fact that both atoms are in the ground state will result 
corresponding to an enhanced symmetry (see 
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Methods). While experimentally more challenging, the use 
of two different ground state species will solve the prob- 
lem of losses associated with collisions of two excited state 
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FIG. 2: Young diagrams describing the irreducible representa- 
tions of SU(N) on individual sites, a, A general diagram consists 
of rij boxes arranged into at most two columns (to satisfy fermionic 
antisymmetry with only two orbital states) whose heights we will 
denote by p and g, such that N > p > q and p + q = rij. See 
Supplementary Information for a brief review of Young diagrams, b, 
The Young diagrams for the two special cases discussed in the main 
text: (1) {p,q) = (1, 0) and (2) (p, g) = (p, 0) on a bipartite lattice. 



atoms and will reduce the (already very weak) nuclear-spin- 
dependence of flee and Ueg. 



Spin Hamiltonians 

One of the simplest interesting limits of Eq. (|2|i is the 
strongly interacting regime (J/U <C 1) where the Hilbert 
space is restricted to a given energy manifold of the Jg = 
Je = Hamiltonian (with a fixed number of atoms on each 
site), and tunneling is allowed only virtually, giving rise to 
an effective spin (and pseudo-spin) Hamiltonian. Single-site 
energy manifolds can be classified according to the number 
of atoms rij = rijg + nje, the pseudo-spin component Tj, 
and the spin symmetry (SU(N) representation) described by 
a Young diagram. As shown in Fig. |2^, each diagram con- 
sists of rij boxes and at most two columns of heights p and q, 
representing two sets of antisymmetrized indices. 

The [/(I) X SU{N) symmetry of Eq. (|2]i restricts the order 
spin Hamiltonian to the form 



(id), a 

+ E yen^gn,e + X%Sl{i, 9)Sl\j, e) 



where the sum over n and m is implied in all but the n terms 

and 5*^,70) = '^]i3n^jam- {i ^ i] means that all 4 preceding 
terms are repeated with i and j exchanged. The coefficients k, 
A, K, and A are of order ,P /U with the exact form determined 
by what single-site energy manifolds we are considering, k 
terms describe nearest neighbor repulsion or attraction, while 
A, K, and A terms describe nearest neighbor exchange of spins, 
pseudo-spins, and complete atomic states, respectively. With- 
out loss of generality, kJ^ = k^' and A^ = A^'. In many 
cases (e.g. case (2) below), the Hilbert space, which i?(p,<j) 
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FIG. 3: The ground-state phase diagram for the SU(N=2) Kugel- 
Khomskii model restricted to two wells, left (L) and right (R). a, 

The phase diagram for T~ — —1 (two g atoms). \gg) = \gg)LR- 
\s) and \t) are spin singlet and triplet states, respectively, b, The 
phase diagram for T, = (one g atom and one e atom). |E) = 
75(|eff>i« - l3e>ifl) and |r) = -^{\eg)LR + \ge)LR) are anti- 
symmetric and symmetric orbital states, respectively. See Supple- 
mentary Information for a detailed discussion of both of these dia- 
grams. 



acts on, has riie and riig constant for all i, which not only 
forces Kg^^ = Ag-'g = but also allows one to ignore the con- 
stant K^J and terms. We now discuss two special cases of 
^{p,g} shown in Fig.|2]3. A third case, (p, q) ~ (1, 1), which 
reduces for = 2 to the spin-1 Heisenberg antiferromagnet 
is discussed in the Supplementary Information. 

(1) In the case of one atom per site, {p,q) = (1,0). H(^p g-j is 
then a generalization to arbitrary N of the SU(A^ = 2) Kugel- 
Khomskii model lUdll], and we rewrite it as (see Supplemen- 
tary Information) 



-[A + BSf^](TtT^ + -) + h{l - Sf^)(Tt + 7- 



, (6) 



where Sfj = S';^(i)^;r(i) is +1 (-1) for a symmetric 

(antisymmetric) spin state, A = 2k 



Xg, B — 2Age 



Ke + Kg, and h = — Hg)/2. The N ^ 2 Kugel-Khomskii 
Hamiltonian is used to model the spin-orbital interactions (not 
to be confused with relativistic spin-orbit coupling) in transi- 
tion metal oxides with perovskite structure fl3\. Our imple- 
mentation allows to realize clean spin-orbital interactions un- 
altered by lattice and Jahn-Teller distortions present in solids 

To get a sense of the competing spin and orbital orders 
ifiol [III [12I1 characterizing i/(i.o), we consider the simplest 
case of only two sites {L and R) and N = 2 (with spin states 
denoted by j and l). To avoid losses in e-e collisions, we set 
Uee = 00 (see Supplementary Information). The double-well 
ground-state phase diagram for = 1 (two e atoms) is then 
trivial, while the = — 1 (two g atoms) and = (one g 
atom and one e atom) diagrams are shown in Fig. [3] One can 
see that, depending on the signs and relative magnitudes of the 
interactions, various combinations of ferromagnetic (triplet) 
and antiferromagnetic (singlet) spin and orbital orders are fa- 
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FIG. 4: Probing the phases of the SU(N) antiferromagnet on a 2D 
square lattice, a shows the phase diagram for the case ha + riB = 
A^. Some points on this diagram have been explored in earlier nu- 
merical studies t2^[30l . [3Tll and are marked according to the ground 
state obtained: Neel (circles), columnar-valence-bond solid (VBS) 
[shown schematically in b] (squares), and possibly critical spin liq- 
uid (triangle) (30l[3lll . Since for sufficiently large N quantum fluc- 
tuations tend to destabilize long-range magnetic ordering, it is likely 
that VBS ordering characterizes the ground state for all A > 4 (i.e. 
above the wavy line). 



vored. In the Methods, we propose a double-well experiment 
along the lines of Ref. ll39ll to probe the spin-orbital interac- 
tions giving rise to the = diagram in Fig.[3j5. Multi-well 
extensions of this experiment may shed light on the model's 
many -body phase diagram, which has been studied for N = 2 
and mostly at mean-field level or in special cases, such as in 
one dimension or in the presence of enhanced symmetries (see 

e.g. iHlnllll). 

(2) In order to study SU(N) spin physics alone, we con- 
sider the case of g atoms only. On a bipartite lattice with 
sublattices A and B, we choose A sites to have ua < N 
atoms [{p, q) — {ua, 0)] and B sites to have ns < N atoms 
[(P: 9) = ('T-BiO)]. This setup can be engineered in cold atoms 
by using a superlattice to adjust the depths of the two sublat- 
tices favoring a higher filling factor in deeper wells. Hf^p g) 
then reduces to 



H{p,0) - 



2 T^U 
■^'^g'-^gg 



^99 



[UggiuA - Ub) + Ay 



(7) 



where A is the energy offset between adjacent lattice sites. 
The coupling constant can be made either positive (antifer- 
romagnetic) or negative (ferromagnetic) depending on the 
choice of parameters |39]. Three body recombination pro- 
cesses will likely limit the lifetime of the atoms when nj > 3 
(see Supplementary Information). 

We focus on the 2D square lattice in the antiferromag- 
netic regime. The case ua + ns ~ N shares with the 
SU(2) Heisenberg model the crucial property that two ad- 
jacent spins can form an SU(N) singlet, and has thus been 
studied extensively as a large-N generalization of SU(2) mag- 
netism ll27i I28I1 . Fig. shows the expected phase diagram 
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for the case ua + nB = N, which features Neel (circles), 
valence-bond-soHd (VBS) (s quar es) [Fig. [^], and possible 
critical spin liquid (triangle) OOi 13111 ground states. To ac- 
cess various ground states of the system, the initial state must 
be carefully prepared so that the conserved quantities 5™ take 
values appropriate for these ground states. Another interesting 
and experimentally relevant case, ha = ns N/2, which 
can also exhibit spin liquid and VBS-type ground states, is 
discussed in the Supplementary Information and in Ref. Isill . 

Since one can vary N just by choosing the number of ini- 
tially populated Zeeman levels (e.g. via a combination of opti- 
cal pumping and coherent manipulation), alkaline-earth atoms 
offer a unique arena to probe the phase diagram of i/(p.o)- in- 
cluding exotic phases such as VBS [Fig.|33], as well as com- 
peting magnetically ordered states. We propose to load a band 
insulator of N g atoms per site, then slowly split each well into 
two to form an array of independent SU(N) singlets in a pat- 
tern shown in Fig. |4}5. The intersinglet tunneling rate should 
then be adiabatically increased up to the intrasinglet tunneling 
rate. As N increases, the magnetic or singlet nature of the 
state can be probed by measuring the Neel order parameter 
(see the description of the Kugel-Khomskii double-well ex- 
periment in the Methods) and spin-spin correlations via noise 
spectroscopy in the time of flight 140(1 (which directly mea- 
sures E,, 5)^;;(J, <7))e^«(*-^)). 



The Kondo lattice model (KLM) 

The SU(N) Kondo lattice model (KLM) |[il[T3l is another 
example of the rich physics, beyond the Mott regime, which 
could be simulated with alkaline-earth atoms. The KLM is 
one of the canonical models used to study strongly corre- 
lated electron systems, such as manganese oxide perovskites 
ll20ll and rare earth and actinide compounds classed as heavy 
fermion materials ll25ll . 

For its implementation with cold atoms (for iV = 2, see 
also Refs. ll23il24ll ). we propose to put one e atom (localized 
spin) per site in a deep lattice such that Je <C Uee, so that 
we can set J,. = and rije ~ 1 for all j in Eq. (|2]). We also 
suppose that we can set Ugg ~ 0, e.g. by taking a very shallow 
q-lattice (see Fig.|5^). The resulting Hamiltonian is the SU(N) 
KLM dill 




^ex \l J g 



HKLM--^Jg{c\^ 
{j,i)m 



The magnitude of Vex can be adjusted by shifting the e and g 
lattices relative to each other |01 ■ 

The properties of the SU(N) KLM depend crucially on the 
sign of the exchange interaction. For concreteness, we focus 
on the antiferromagnetic (AF) case {Vex < 0), which favors 
formation of spin-antisymmetric states (singlets, for iV = 2) 
between mobile fermions and localized spins. This regime 
describes the physics of heavy fermion materials ll25ll . and, 
in the case of a single localized spin, gives rise to the Kondo 
effect. 

In the limit | Vex \ Jg, 9 atoms mediate long-range RKKY 



FIG. 5: Kondo lattice model for the case N = 2. a. The schematic 
of the setup, g atoms are green; e atoms are yellow; the spin ba- 
sis is {t, J,}, b, Schematic representation of the competition be- 
tween RKKY magnetism vs Kondo singlet formation in the SU(2) 
AF KLM (see Jig . |2^ and references therein). In this model, 
the localized spin- 1/2 e atoms couple antiferromagnetically to the 
delocalized g atoms, via an on-site exchange interaction propor- 
tional to Vex ■ This coupling favors the formation of localized Kondo 
singlets between e and g atoms, with characteristic energy scale 
ksTx ~ Jg exp(— cJt,/|V"e3:|), with c a dimensionless constant of 
order one 1231 . On the other hand, the g atoms can mediate long- 
range RKKY interactions between the e atoms, giving rise to mag- 
netic order (which can be antiferromagnetic (AF) or ferromagnetic 
depending on the density of g atoms), where the characteristic en- 
ergy is ksTRKKY ~ Vex/ Jg- The competition between Kondo ef- 
fect and RKKY magnetism leads to very rich physics. For small val- 
ues of \Vex\/ Jg, the RKKY interaction is dominant and the system 
orders magnetically. At intermediate values of \Vex\/ Jg, the energy 
scales Tk and Trkky are of comparable strength, and a variety of 
novel quantum phenomena are expected to arise, including quantum 
criticality and non-Fermi liquid (NFL) physics ll25ll2^ . With further 
increase of the \Vex\/ Jg coupling, magnetic order is suppressed, the 
localized e atoms become screened into singlet states and melt into 
>-)the g-atom Fermi sea, forming the so called heavy Fermi liquid state 
(HFL). The large Fermi volume I21I1 . which is set by the total num- 
ber of g atoms plus e atoms, can be directly probed by measuring the 
momentum distribution via time of flight imaging. 



interactions 111411 between localized spins and tend to induce 
magnetic ordering (antiferromagnetic or ferromagnetic de- 
pending on the density of g atoms) of the latter, at least for 
N ^ 2. The engineering of RKKY interactions can be tested 
in an array of isolated double wells (see Methods). At in- 
termediate and large \Vex\, the formation of Kondo singlets 
dominates the RKKY interaction and favors a magnetically 
disordered heavy Fermi liquid (HFL) ground state with a sig- 



6 



nificantly enhanced effective quasiparticle mass (see Fig.lSj)). 
The competition between RKKY interactions and the Kondo 
effect in the regime where both are comparable is subtle, and 
the resulting phases and phase transitions ll25i l26ll are not 
well-understood. Ultracold alkaline-earth atoms provide a 
promising platform to study these phases and phase transi- 
tions. 

In the large-iV Umit lillil, the SU(N) HFL can be con- 
trollably studied, and 1/N expansions have successfully re- 
produced the experimentally observed properties of the HFL. 
However, very little is known about the SU(N) model outside 
the HFL regime. Several very interesting parameter regimes 
in this domain can be directly probed with our system, as dis- 
cussed in the Methods. 



Experimental Accessibility 

The phenomena described in this manuscript can be probed 
with experimental systems under development. Indeed, we 
show in the Methods that SU(N)-breaking terms are suffi- 
ciently weak, and here we discuss the temperature require- 
ments. 

The key energy scale in the spin Hamiltonians [Eq. ©I 
is the superexchange energy J^/[/, while the RKKY en- 
ergy scale is ksTjiKKY ~ ^ex/Jg- I" their region of va- 
lidity (J < U and \Vex\ < Jg, respectively), these en- 
ergy scales are limited from above by the interaction energy 
{U and |Ve2-|, respectively), which typically corresponds to 
temperatures T < lOOnK ll39ll . Thanks to the additional 
cooling associated with certain adiabatic changes ll4ll l42ll . 
T ^ lOnK and the Mott insulating regime have already been 
achieved with fermionic alkali atoms ll43ll . and are therefore 
expected to be achievable with fermionic alkaline-earths, as 
well (a bosonic alkaline-earth Mott insulator has already been 
achieved ||44l] ). Furthermore, the requirement to reach ksT 
smaller than J^/[/ or V^^/Jg can often be relaxed. First, 
the double-well experiments, such as the ones discussed in 
the Methods in the contexts of the Kugel-Khomskii and the 
Kondo lattice models, are performed out of thermal equilib- 
rium, and can, thus, access energy scales far below the tem- 
perature of the original cloud |39i|. Second, for SU(N) anti- 
ferromagnets, the energy range between J^/[/ and N /U 
may also exhibit intriguing physics: in this regime, SU(N) 
singlets, which require N,P /U energy to break, stay intact 
but can diffuse around. Finally, in the Vex < Kondo lattice 
model, exotic heavy Fermi liquid behavior is expected when 
Jg < I Vex I and the temperature is below the Kondo tempera- 
ture, i.e. fcsT < Jg cxp(— cJg/|Ve.T|) with c is a dimension- 
less constant of order one |25i] . Thus, with Jg chosen to be on 
the order of | Vex | , fcs T as high as ^ \Vex \ may be sufficient. 



Outlooli 

The proposed experiments should be regarded as bridges 
aiming to connect well-understood physics to the complex and 
poorly understood behavior of strongly correlated systems. It 



is important to emphasize that, except for the one dimensional 
case, the phase diagram of most of the models considered is 
only known at mean field level or numerically in reduced sys- 
tem sizes. Therefore, their experimental realization in clean 
and controllable ultracold atomic systems can provide major 
advances. 

Our proposal motivates other new lines of research. Ultra- 
cold bosonic or fermionic diatomic molecules |45] may give 
rise to similar SU(N) models with large N and with the possi- 
bility of long-range interactions. Ions with alkaline-earth-like 
structure, such as A1+ could also be considered in this context. 
It would also be interesting to explore the possibility of real- 
izing topological phases with SU(N) models for applications 
in topological quantum computation ll34ll . Beyond quantum 
magnetism, the fact that the formation of SU(N) singlets re- 
quires N partners might give rise to novel exotic types of su- 
perfluidity and novel types of BCS-BEC crossover IBTIl . Prac- 
tical applications of our Hubbard model, such as the calcula- 
tion of the collisional frequency shift in atomic clocks ll46ll . 
can also be foreseen. 



Note added in proof. After the submission of this paper, a 
theoretical study of the SU(6)-symmetric ^^'^Yb system was 
reported ifsoll . 



Metliods 

Experimental tools available for alkaline-earth atoms 

Many experimental tools, such as tuning the interaction 
strength by adjusting laser intensities ll39ll . are common to 
both alkali and alkaline-earth atoms. There are, however, 
some experimental tools specific to alkaline earths; we review 
them in this Section. 

First, a combination of optical pumping 101 and direct co- 
herent manipulation of the \g) — |e) transition in the presence 
of a magnetic field ([ll01 can be used IH] to prepare any desired 
single-atom state within the 2 (21+ l)-dimensional manifold 
with basis \am), where a = g ot e and rn = — /, ...,/. This 
coherent manipulation can also be used to exchange quantum 
information between nuclear spin states and electronic states. 
Second, by using far-detuned probe light or a large magnetic 
field to decouple the electronic angular momentum J and the 
nuclear spin /, the electronic \g) — |e) degree of freedom can 
be measured by collecting fluorescence without destroying the 
nuclear spin state [Isl]. Fluorescence measurement of the nu- 
clear spins can be achieved by mapping nuclear spin states 
onto electronic states ItIHI: for example, for a spin- 1/2 nu- 
cleus, a TT pulse between \g,m ~ 1/2) and |e,?T!, = —1/2) 
allows one to accomplish a swap gate between the nuclear 
{1/2, —1/2} qubit and the electronic {e, g} qubit. Single-site 
spatial resolution during the coherent manipulation and fluo- 
rescence measurement can be achieved using magnetic field 
gradients |01 or dark-state-based techniques lH H?!] that rely 
on an auxiliary laser field whose intensity vanishes at certain 
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locations. Third, an appropriate choice of laser frequencies 
allows one to obtain independent lattices for states g and e 
JtII- Finally, optical Feshbach resonances ll48ll may be used to 
control scattering lengths site-specifically and nearly instanta- 
neously. 



only and extending labels m and n to run over 2N states in- 
stead of N . 



Double-well Kugel-Khomskii and RKKY experiments 



Enhanced Symmetries 

While in the general case, our Hubbard model [Eq. ^] 
satisfies U{1) x SU{N) symmetry, for particular choices of 
parameters, higher symmetry is possible. In particular, if 
Jg = Je and the interaction energies for all states within the 



pseudo-spin triplet are equal (U, 
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U+), the full 



SU(2) symmetry (not just U(l)) in the pseudo-spin space is 
satisfied. Alternatively, if Vex ~ 0, then both S™{i,g) and 
S™{i, e) generate SU(N) symmetries resulting in the overall 
U{1) X SU{N) X SU{N) symmetry. Finally, if both condi- 
tions are satisfied, i.e. all four Ux are equal and Jg — J^, then 
H satisfies the full SU(2N) symmetry (2N can be as high as 
20) generated by 



nam 
^f3n 



j 

U ■ 



in which case the interaction reduces to ^ V. 



(9) 



1), 



where Uj = rijg 



In the case when |e) and \g) correspond to two ground states 
of two different atoms (with nuclear spin /g and Ig, respec- 
tively), we will have = a~g (i.e Vex — 0), which is equiva- 
lent to imposing t/(l) x S';7(Afg = 2Ig + l)-KSU{Ne = 2Ie + 
1) symmetry, where SU{2Ia + 1) is generated by S'^{i, a). 
While for Ig ^ le, the m index in Cjam will run over a dif- 
ferent set of values depending on a, the Hubbard Hamiltonian 
will still have the form of Eq. (|2|i (except with Vex = 0). If 
one further assumes that Jg = Je and Ugg = Uee = Ueg, 
the interaction satisfies the full SU{Ng + Ne) symmetry. It 
is worth noting that for the case of two different ground state 
atoms, this higher symmetry is easier to achieve than for the 
case of two internal states of the same atom, since a+ = a^g 
automatically. Thus, in particular, it might be possible to ob- 
tain 5f/(18) with 87Sr (/ = 9/2) and ^^Ca (I = 7/2) simply 
by adjusting the intensities of the two lattices (to set Jg = Je 
and Ugg = Uee) and then shifting the two lattices relative to 
each other (to set Ueg = Ugg). 

Enhanced symmetries of the Hubbard model [Eq. (|2]l] are 
inherited by the spin Hamiltonian [Eq. (|5]l]. In particular, im- 
posing SU{2) X SU{N) instead of C/(l) x SUiN) forces 



Ne 



ge ge' 



X^3 = X^J = X^3 



\^e = A^J. 



Alternatively, imposing U{1) x SU(N) x SU(N) forces 



"ge 



= Ag-'g = 0. Finally, imposing the full SU(2N) forces 



the satisfaction of both sets of conditions, yielding 



(10) 



which is, of course, equivalent to restricting Eq. (|5]l to g-atoms 



In the main text and in the following Methods Section, we 
discuss the open questions and previously unexplored regimes 
associated with the SU(N) Kugel-Khomskii and Kondo lattice 
models (KLM) that can be studied with ultracold aUcaline- 
earth atoms. As a stepping stone toward these many-body ex- 
periments, we propose in this Section two proof-of-principle 
experiments in an array of isolated double wells with N ^ 2 
(with the spin basis {f, j}): one to probe the spin-orbital in- 
teractions of the Kugel-Khomskii model and one to probe the 
RKKY interactions associated with KLM. 

We first propose an experiment along the lines of Ref. ll39ll 
to probe the spin-orbital interactions giving rise to the = 
diagram in Fig.|3j). In the Supplementary Information, we de- 
scribe how to prepare an array of independent double wells in 
the state |e, T)-l|.9i i)-R' which is a superposition of the four 
eigenstates featured in Fig. |3j)- The energies of these four 
eigenstates [Eqs. ( IS4HS7l i1 can be extracted from the Fourier 
analysis of the population imbalance as a function of time: 



AN{t) = 



neR 



T-gR — neL 



hu^c 



COS 



AN can be measured by combining the dump- 
ing technique, band mapping, and Stern-Gerlach filtering of 
Ref. ll39ll with the use of two probe laser frequencies to distin- 
guish between \g) and |e). 

We now turn to the double-well experiment aimed at prob- 
ing RKKY interactions. After preparing the state -^[\g,l 

)l + Is, i)j?)|e, i)L|e, t)i?, (see Supplementary Information 
for how to prepare this state), we propose to monitor the 
Neel order parameter for the e atoms, Nez = \[ne]L — 



nelL 
1 



COS 



{ne]R 



UelR)]. In the limit \Vex\ < Jg, Nez{t) 



3 — V^) - I COS - 1^ j [in the Supplemen- 

tary Information, we present the plot of Nez{t) for Vex = 
— Jg/10]. It exhibits fast oscillations with frequency ~ Vex, 
modulated by an envelope of frequency ~ V^^ j Jg induced 
by RKKY interactions. In order to probe RKKY interactions 
only, it is important to suppress super-exchange ^Je/Uee and 
thus to choose Je/Uee Small. To study the full spatial depen- 
dence of RKKY interactions, one must of course go beyond 
the double-well setup. We also note that recent experiments 
using alkali atoms populating the lowest two vibrational levels 
of a deep optical lattice have measured the local singlet-triplet 
splitting induced by Vex ll49ll . 



3 Kit 



Physics accessible with the alkaline-earth Kondo lattice model 

The alkaline-earth atom realization of the AF KLM is well- 
suited to access a number of parameter regimes that are out 
of reach in solid state materials. One example is the one 
dimensional (ID) limit, since, to our knowledge, real solid 
state materials exhibiting KLM physics are restricted to 2D or 
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3D. Another example is the regime of large Kondo exchange 
{\Vex\ ^ Jg)^ which is interesting even for N = 2. In this 
limit the system is well described by the J7 — > oo Hubbard 
model ifisll by identifying the Kondo singlets with empty sites 
(holes) and the unpaired localized spins with hard core elec- 
trons. From this mapping, possible ferromagnetic ordering is 
expected at small hole concentration (small rig), however the 
stability of this phase for increasing hole concentration and 
finite \Vex\ values remains unknown. For general N, in the 
extreme limit Jg = 0, the ground state is highly degenerate: 
for any distribution of the g atom density njg < N, there is 
a ground state (with further spin degeneracy), where on each 
site the spins combine antisymmetrically to minimize the ex- 
change interaction. Lifting of such extensive degeneracies of- 
ten leads to novel ground states; this will be addressed in fu- 
ture studies using degenerate perturbation theory in Jg/Vex- 
For N > 2, AF SU(N) spin models have a different kind of 
extensive degeneracy, which was argued to destroy antiferro- 
magnetism and to lead to non-magnetic spin liquid and VBS- 
like ground states ll34ll . Similar expectations are likely to ap- 
ply to the KLM at small iVexl/Jg, where the = 2 antifer- 
romagnetism may give way to situations where the localized 
spins form a non-magnetic state that is effectively decoupled 
from the mobile fermions ll22ll . 

Even though we have set Ugg to zero in Eq. dHJ, it can be 
tuned, for example, by adjusting the g-lattice depth and can 
give rise to interesting physics. For example, the Ug = 1 
case, which is known to be for N = 2 either an antifer- 
romagnetic insulator or a Kondo insulator depending on the 
ratio I Vex I /"/g !19], will become for large enough Ugg and 
A^ > 2 a Mott insulator, because the two atoms on each site 
cannot combine to form an SU(N) singlet. If Ug is reduced 
from unity, the doping of this Mott insulator can be studied, 
and it will be interesting to understand how this physics, usu- 
ally associated with cuprate superconductors, is related to the 
other ground states of the KLM, usually associated with heavy 
fermion compounds. 

Experimental Accessibility 

Immediate experimental accessibility makes our proposal 
particularly appealing. Having shown in the main text that the 



temperature requirements of our proposal are within reach of 
current experimental systems, here we show that the nuclear- 
spin dependence of interaction energies is sufficiently weak to 
keep the SU(N) physics intact. 

In the Supplementary Information, nuclear-spin-dependent 
variation in the interaction energies is estimated to be 

AUgg/Ugg ^ 10'^ HUd AUee/Uee ^ AU^g/U^g ^ 10-3. 

Since the scale of SU(N) breaking is at most AU, a very con- 
servative condition for the physics to be unaffected by SU(N) 
breaking is that all important energy scales are greater than 
AU. In particular, in the spin models with more than one atom 
per site, the condition is AU <^ J'^/U, which can be satisfied 
simultaneously with J <^ U even for AU /U ~ IQ-'^. With 
one atom per site, the SU(N) breaking scale is not AU but 
rather ( J/C/)^A[/, which relaxes the condition to the imme- 
diately satisfied AU /U <C 1. Similarly, in the Kondo lattice 
model, the conditions AVex ^ >/, \ Vex\ can be satisfied for 
AV,^/\Vex\ - 10-3. 
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SUPPLEMENTARY ONLINE MATERMLS 

Known Scattering Lengths 

Very few scattering lengths ax {X = gg, ee, eg"*", eg~) be- 
tween g CSq) and e CPq) states of aUcaUne-earth-like atoms 
are known at the moment, ag n is known for all isotopic com- 
binations of Yb ifsill and Sr il52ll . Estimates of a^e for ®^Sr 
ifssll and of a^g for ®^Sr i54ll also exist. Finally, there is a 
proposal describing how to measure a+g via clock shifts ifssll . 



Nuclear-Spin Independence of the Scattering Lengths 

Independence of scattering lengths from the nuclear spin is 
a key assumption of the paper. This feature allows us to ob- 
tain SU(N)-symmetric models with N as large as 10 and dis- 
tinguishes alkaline-earth atoms from alkali atoms, which can 
exhibit at most an SO(5) symmetry ll56ll57il58ll59ll . a symme- 
try that is weaker than SU(4). The assumption of nuclear-spin 
independence of scattering lengths is consistent with recent 
experiments, where - within experimental precision - the clock 
shift does not depend on how the Zeeman levels are populated 
ll60l l6lll . In this Section, we present the theoretical justifica- 
tion of this assumption. 

Direct magnetic dipole-dipole coupling between the nu- 
clear spins of two atoms sitting on the same site of an op- 
tical lattice is negligible: even for two magnetic dipole mo- 
ments as large as 10 nuclear magnetons at a distance of 10 nm 
(which is significantly smaller than the confinement typically 
achieved in optical lattices ll62ll ). the interaction energy still 
corresponds to a frequency smaller than one Hertz. Therefore, 
nuclei can affect the collisions only via the electrons. All four 
scattering lengths (ogg, afg, and flee) are, thus, expected to 
be independent of the nuclear spin because both g and e have 
total electronic angular momentum J equal to zero, which re- 
sults in the decoupling between nuclear and electronic degrees 
of freedom during the course of a collision. The decoupling 
during a collision is a consequence of the fact that each of the 
four molecular electronic states that correlate with the J = 
separated atom pair has zero projection of total electronic 
angular momentum on the molecular axis. The nuclear spins 
in this case can only couple very weakly to other molecular 
states, even if there is a molecular curve crossing. 

While the short-range potential energy structure for a 
molecule like Sr2 is very complex for the excited states 
ll63il64ll . we will now show that scattering length differences 
among different combinations of nuclear spin projections for 
the same = potential are expected to be very small. The 
scattering length a can be computed as a = a[l — tan($ — 
7r/8)], where a is the average scattering length governed by 
the asymptotic behavior of the potential and $ is the semiclas- 
sical phase computed at zero energy from the classical turn- 
ing point Ro to infinity: $ dRy/ M[~V{R)]/h, where 
—V{R) is the (positive) depth of the interaction potential at 
separation R and M/2 is the reduced mass ll65ll . Defining 
R{t) as the classical trajectory from time t = — oo to time t = 



oo of a particle of mass M/2 at zero energy in the potential 
V{R), we can rewrite the phase as $ = — dtV{R{t))/ h. 
The order of magnitude of the change in the phase asso- 
ciated with different nuclear spin projections can, thus, be es- 
timated as (5$ ~ AtSV/h, where At is the total time in the 
short-range part of the collision and 5V is the typical energy 
difference associated with different nuclear spin projections 
during this time. Since 5V vanishes at i? ^ oo, only the short 
range molecular region contributes to the phase difference. 
Therefore, assuming 5$ ^ 1, a ^ a, and | cos(<I> — 7r/8)| ^ 
1, the nuclear-spin-dependent variation Sa in the scattering 
length can be estimated as Sa/a ^ ^ AtSV/h. 

Turning to the actual numbers. At can be estimated from 
the depth (^ lO'^cm" ^/ic) and the range 10 Bohr radii) 
of the appropriate interatomic potential (see e.g. ll63i l64ll ) to 
be At w 1 ps. For g-g collisions, 5V/h can be estimated 
by the second-order formula E^^/ (hEopi) ^ 200 Hz, where 
Ehf/h ^ 300MHz is the approximate value for the hyper- 
fine spUttings in ^Pi in ^''Sr and Eopt/h 400 THz is the 
optical energy difference between ^Sq and "^Pi in ®^Sr. This 
yields the following estimate for the dependence of agg on the 
nuclear spin: Sugg/agg ^ 5^ ^ 10^^. For e-e and e-g colli- 
sions, an analogous second-order formula would use the fine 
structure splitting between "^Pi and '^Pq in ^^Sr (Ef/h ~ 6 
THz) instead of i^opt to yield ^ 10^*^. However, the latter 
estimate ((5$ ~ 10~^) is too optimistic since molecular states 
that are split by Ef at large interatomic separations may come 
orders of magnitude closer at short range ll66ll . Therefore, a 
more realistic conservative estimate would use the first-order 
formula SV ~ E^f to yield Saee/a-ee ^ ^ 5'^ ^ 

It is important to note, however, that these are all only 
very rough estimates. For example, hyperfine coupling in a 
molecule will differ from the hyperfine coupling in separated 
atoms. In fact, since it is very difficult to predict 5a/ a accu- 
rately, these values would need to be measured. To conclude 
this Section, we would like to emphasize that, as mentioned 
in the main text, if the small nuclear-spin dependence of a^e 
and a^g is not negligible for some applications, one can use 
two different ground state atomic species instead of a ground 
and an excited state of one species. 



Likelihood of Lossy e-e Collisions and Possible Solutions 

Collisions of two e atoms are likely to be accompanied by 
large loss ifssll . This can occur if the molecular 0+ potential 
that correlates with the e-e atoms undergoes an avoided cross- 
ing with a potential curve that correlates with a lower energy 
pair of separated atoms (see, for example, Ref. ll64ll ). Simi- 
lar crossings that result in inelastic energy transfer collisions 
were examined for ^Pi+^Sq collisions of alkaline earth atoms 
in Ref. ||67|l. The likelihood of a relatively high probability of 
an inelastic event during such a crossing with species such 
as Sr or Yb means that the imaginary part hee of the scatter- 
ing length is expected to be large. However, just like aee, b^e 
can not be calculated accurately from the potentials but would 
need to be measured. 

The possible effects of bee on the four examples we discuss 
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(p,9) = (l,l) 



FIG. S2: (p,q) = (1,1) Young diagram. 



FIG. SI: A general Young diagram. 



[Eqs. ( I6ll8] ) and Eq. dsB l are as follows. i?(p,o) [Eq. ©] is, 
of course, not affected because it involves only g atoms. In 
1) [Eq. ( ISlb 1 and Hklm [Eq. ©l, the e lattice is assumed 
to be so deep that is negligible compared to Uee + Vex and 
t/ee, respectively, or to the experimental timescale, thus, fully 
suppressing tunneling of e atoms and occupation of one site 
by more than one e atom. The presence of an imaginary part 
hee of the e-e scattering length will give an effective nonzero 
width to the state with more than one e atom per site and can, 
therefore, only further suppress this tunneling by a Zeno-like 
effect iHm-IMl- 

Therefore, ^f(i,o) [Eq. Q] is the only example that can be 
affected by large hee- In order for i?(i,o) to contain a non- 
negligible term proportional to J? /Uee, the ratio \bee/aee\ 
would need to be very small ItiIi . Several approaches to 
avoid the losses associated with bee in ^(i.o) ^re possible. 
First, the large variety of stable atoms with two valence elec- 
trons (which includes not only alkaline-earths, but also Zn, 
Cd, Hg, and Yb) may have coincidentally an isotope with 
small \bee/aee\, which is more likely for lighter atoms ll67ll . 
Second, while obtaining a good optical Feshbach resonance 
IUIItIIItIItUtI to reduce \bee/aee\ might not be possible, 
it should be possible to use optical Feshbach resonances to 
enhance bee and, thus, suppress ll68l l69i [toIi the virtual occu- 
pation of one site by two e atoms; i?(i,o) would then have the 
same form as in Eq. (|6]l, except with Uee effectively set to in- 
finity. Notice that here we suggest to use optical Feshbach 
resonances to affect e-e scattering, which is different from 
the typical apphcation to g-g scattering ll5lll72ll73ll74ll75ll . 
Third, one can consider using a different ground state atom to 
represent state |e), which would set Vex = in if(i,o)- Fi- 
nally, one could simply use an e-lattice that is deep enough to 
make Je negligible, which would, however, lead to the loss of 
terms in i?(i.o) that exchange the pseudospin between neigh- 
boring sites. 



Brief Review of Young Diagrams 

Irreducible representations of SU(2) are classified accord- 
ing to the total half-integer angular momentum J and have di- 
mension 2 J + 1. On the other hand, a (semistandard) Young 
diagram, instead of a single value J, is used to describe an 
irreducible representation of SU(N) for a general N ll76l[77ll . 
As shown in the example in Fig. SI, a Young diagram has all 
its rows left-aligned, has the length of rows weakly decreasing 
from top to bottom, and has at most rows. The dimension 
of the representation corresponding to a given diagram is the 



number of ways to fill the diagram with integers from 1 to A^ 
such that the numbers weakly increase across each row and 
strictly increase down each column. For our purposes, the 
number of boxes in the diagram is the number of atoms on the 
site, and the diagram describes the (nuclear) spin symmetry of 
the particular chosen single-site energy manifold. In particu- 
lar, columns represent antisymmetrized indices, while rows 
are related to (but do not directly represent) symmetrized in- 
dices. It is the relation between antisymmetrized indices and 
the columns that limits the number of rows to A. On the other 
hand, since the full wavefunction (spin and orbital) on each 
site must satisfy complete fermionic antisymmetry, the rela- 
tion between rows and symmetrized indices and the fact that 
we have only two orbital states (g and e) force all our diagrams 
to have at most two columns. 



The (p, q) = (1, 1) spin Hamiltonian and the spin-1 Heisenberg 
antiferromagnet 

In the main text, we discussed two special cases of the 
spin Hamiltonian ff(p.q), both of which had a single-column 
SU(N) representation on each site (i.e. q = 0). In this Sec- 
tion, we discuss the simplest SU(N) representation with two 
columns, {p,q) — (1,1) [see Fig. S2]. It can be obtained 
when there is one g and one e atom per site in the electronic 
singlet \ge) — \eg) configuration. Setting Jg = to avoid e-e 
collisions, reduces to 



(1.1) 



(SI) 



The case of A^ = 2 is the spin-1 antiferromagnetic Heisen- 
berg model. This model has a ID ground state with hidden 
topological structure ItsIi . Recently, applications of related 
models in one-way quantum computation have been proposed 
||79[ Isoll . Models with more complicated two-column repre- 
sentations may have exotic chiral spin liquid ground states that 
support non-Abelian anyons and that might thus be used for 
topological quantum computation llsill . 



The Kugel-Khomskii model and the double-well phase diagram 

In the main text, we omitted the values of the parame- 
ters in 7J(p^that characterize the Kugel-Khomskii model 
0) [Eq. ©I. In this Section, we present these parame- 
ters. We also present a detailed discussion of the double-well 
case phase diagram. 



12 



The parameters in Hf^p q-^ that characterize the Kugel 
Khomskii model i^(i.o) [Eq. ©] are A*-' = = = 

-Kg, Ag 



where |S) 



2(7+ 



2J^ ^ 

■ = A„ 



K 



A„ 



2t/+ 

To avoid loss in e-e collisions, 
oo (see 



2c/r, 



we assume for the rest of this Section that [/, 
Supplementary Information for a discussion of losses in e-e 
collisions). 

The nontrivial orbital-orbital, spin-spin, and spin-orbital in- 
teractions in H^i Q-f [Eq. ©J result in competing orders, with 
the actual ground-state order dependent on the parameters of 
the Hamiltonian -ff(i.o)- To get a sense of the possible orders, 
we consider the case iV = 2 (with the spin states denoted by f 
and l) and discuss the double-well problem, with the wells de- 
noted by L (left) and R (right). Due to the large optical energy 
separating e and g, which we have ignored after Eq. ([TJ, the 



three manifolds of constant 



j,{T^ = -1,0,1) 
should each be considered separately. 

The four states in the = 1 manifold, the subspace of 
two e atoms, are |ee)|s) and |ee)|t). Here |ee) = \ee)LR is 
the orbital (or pseudo-spin) state, while \t) ~ \ TDlHi I II 
)lr, ^(I n)LR + I U)lr) and \s) = ti>LK - I it 
)lr.) are the triplet and singlet spin states. Since Uee = oo, 
all four of these states have zero energy and the ground-state 
phase diagram is trivial. 

The four states in the = — 1 manifold (two g atoms) are 
split by i^(i,o) into two energy manifolds: 



\99)\t), E = Q, 
\99)\s 



4J2 

E = 2.. 

Uaa 



(S2) 

(S3) 



Only \gg)\s) can take advantage of the virtual tunneling since 
two g atoms in the triplet spin states cannot sit on the same 
site. Which of the two manifolds is the ground manifold de- 
pends on the sign of Ugg, as shown in the ground-state phase 
diagram in Fig. [3^. It is important to emphasize that for 
Ugg < 0, the subspace of one g atom per site may be sub- 
ject to extra loss down to the lower energy states that have 
both g atoms in the same well. It is also worth noting that the 
diagram is only valid for Jg <C \Ugg\. 

Finally, the eight states in the T~ ~ Q manifold (one g atom 
and one e atom) are split by i?(i,o) into four energy manifolds: 



{Jg + Je? 



eg 



\t)\s), E = - 



|r)|t), E 



\n\s), E=- 



jJg + Jef 

Utg 

{ Jg - Je? 
Ueg 

{Jg-Je? 

Utg 



(S4) 
(S5) 
(S6) 
(S7) 



7f (|eff)Lfl - \ge)LR) and |r) = -^{\eg)LR 
\ge)LR) are anti-symmetric and symmetric orbital states, re- 
spectively. The denominators U~g and U^g in the energies of 
the \t) and |s) states, respectively, reflect the fact that tunnel- 
ing preserves the nuclear spin. At the same time, the ± signs 
in the numerators can be understood by considering the case 
Jg = Je, when all states with overall symmetry under par- 
ticle exchange must have zero energy since for these states 
tunneUng is forbidden due to the Pauli exclusion principle. 
The corresponding ground-state phase diagram as a function 
of the signs and relative magnitude of U'^g and U'^g is shown in 
Fig- Eh- As in the case of the = 1 phase diagram, negative 
interaction energies may lead to increased losses. 



Effects of Three-Body Recombination 

Three-body recombination S [sl HI IM Isl is a pro- 
cess during which three atoms come together to form a di- 
atomic bound state and a single atom, and both final products 
have enough kinetic energy to leave the trap. While in certain 
cases, three-body recombination can be an asset ItoIi . usually 
it results in the loss of atoms and, thus, limits the duration of 
the experiment. For our purposes, we can describe three-body 
recombination by a decay rate 73 iFtoIi resulting in a loss of 
three particles from one site. This rate will likely depend on 
what atomic states are involved and, to the best of our knowl- 
edge, has not yet been measured or calculated for fermionic 
alkaline-earth atoms. 

Out of the four examples [Eqs. (|6][8l) and Eq. ( ISIM that we 
discuss, only H(i,i) [Eq. ( ISIM and U{p,o) [Eq. (|7]i] may be 
affected by three-body recombination [Hklm [Eq. (jSj] as- 
sumes negligible g-g interactions, such as in a very shallow g 
lattice or with a low density of g atoms). In the case of -^(1,1), 
two g atoms and one e atom occupy the same site virtually 
in the intermediate state that gives rise to the second order 
spin Hamiltonian with interaction strength cx Jg/{Ugg + Vex ) ■ 
Thinking of 73 as an effective linewidth for the intermediate 
state, i^(i,i) will be valid and losses small provided that 73 is 
smaller than the effective "detuning" Ugg + Vex- Since scat- 
tering lengths for alkaline-earth atoms lf5lll52ll54ll are compa- 
rable to those for alkali atoms, Ugg + Vex can be on the order 



of several kHz 16211 . At the same time, I/73 for bosonic al- 
kali atoms in deep traps can be on the order of 1 s ll86ll . If 73 
were the same in our case, 73 <C Ugg + Vex would be satis- 
fied. Ways of controlling the interactions via optical Feshbach 
resonances ll5llF72lF73ll74lF75ll may also be envisioned. 

In the case of ^f(p,o) [Eq- O], (n^, Jt-b) = (li 1) does not 
suffer from three-body recombination. {nA^ns) = (1,2) and 
(2, 2) may have three atoms per site virtually. As in the dis- 
cussion of H(i i), provided 73 associated with three g atoms 
per site is smaller than Ugg, these configurations should be 
accessible. For the case {nA,nB) = (1,2), 73 ^ Ugg is 
also acceptable, since it will effectively prohibit the tunnel- 
ing of the atoms to the state with 3 atoms on a site ItoIi . but 
the interaction can still take place through the intermediate 
state, in which an atom from a B site tunnels to an A site 
and back. One can also envision ways to use optical Feshbach 
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FIG. S3: Square lattice valence plaquette solid for = 4. 

When iV = 4 and ua = ~ 1, four sites are required to 
form an SU(4) singlet; these singlets can in turn form the 
schematically shown plaquette-ordered state or a disordered 
phase made of resonant plaquette states iIstIi . 



resonance techniques ifsil [721 FtsIi to induce large 73. To be 
able to resolve the superexchange coupling ^ Jg/ Ugg in cases 
where ua or ub is equal to 3, one must have 73 < Jg/Ugg. 
Given that superexchange coupling can be as high as 1 kHz 
ll62ll . this condition should also be achievable. Although ha 
or ub greater than 3 will result in even shorter lifetimes ll84ll . 
there is a good chance that relatively large ua and n b can be 
achieved: at least, for bosonic alkali atoms in an n = 5 Mott 
insulator state, the lifetime can still be as long as 0.2 s js^. 



The (p, 0) spin Hamiltonian with ua = ns ^ N/2 

In the main text, we focused on one special case of the anti- 
ferromagnetic (p, 0) spin Hamiltonian on a square lattice, that 
with HA + nB = N (where ua and ub denote the number 
of atoms per site on the two sublattices). In this Section, we 
describe another interesting and experimentally relevant case, 
"A = 7^ iV/2 illlllMllillSlalllllll. Poten- 
tial ground states include states built from valence plaquettes 
(Fig. S3) |88b [1^, resonant plaquette states [87], and topo- 
logical spin liquids ll8lll90ll . Valence plaquette states and res- 
onant plaquette states are the natural generalization of VBS 
states and resonant valence bond states (RVB) ||94l] . respec- 
tively; for example, when ua = nB = 1, N lattice sites are 
needed to form a SU(N) singlet. Fig. S3 depicts a square lat- 
tice valence plaquette solid for ua ~ ub ~ I and iV = 4. 
Techniques for detecting some of these phases are discussed 
in Ref. Isill . The experiment described in the main text for 
the case ua + tlb ~ N can also be generalized to probe the 
nA = nB ^ A^/2 phase diagram including exotic phases such 
as valence plaquette solids [Fig. S3], as well as competing 
magnetically ordered states. The main difference is that after 
preparing a band insulator of N g atoms per site, each site 
should be split not necessarily into two sites but into the num- 
ber of sites that is appropriate for the case being considered 
(e.g. 4 for the case shown in Fig. S3). 




L R L R 

FIG. S4: A schematic diagram describing the preparation 
of the double-well state |e, T)L|ff, Dr- 



Double-well Kugel-Khomskii and RKKY experiments 

In the Methods, we have omitted the description of how 
to prepare the initial states for the proof-of -principle double- 
well Kugel-Khomskii and RKKY experiments. We present 
this description in this Section. 

We first describe how to prepare an array of independent 
double wells in the state |e, T)l|<7, 1)r, which we use for the 
proof-of -principle experiment to probe the spin-orbital inter- 
actio ns in the Kugel-Khomskii model 

llOlll . After loading a band insulator of \g, |) atoms in a deep 
optical lattice, an additional lattice for both g (green) and e 
(yellow) atoms with twice the spacing of the first lattice is 
turned on in one direction to create an array of independent 
double wells ll62ll . Then, as shown in Fig. S4, in the presence 
of an e-lattice bias, (T+ polarized light on resonance with the 
\g, 1)l \e, transition can be used to prepare the state 
|e, T)l|5, i)_R- For examples of earlier orbital physics stud- 
ies with ultracold atoms, where the orbitals are distinguished 
only by the different motional states of the ato ms, we refer the 
reader to Refs. 

iM [M [M [M [M [io3] 

and references 

therein. 

We now describe how to prepare the initial state -^{\g, i 

)l + 1.9, i)i?)|e, J,)L|e, T>i?. (see Fig. S5a) for the double- 
well proof-of-principle RKKY experiment, whose expected 
Neel order parameter Nez{t) for Vex = — Jg/10 is show in 
Fig. S5b. The first step to prepare the initial state -ijdg, | 

)l + 1.9, i)fi)|e, DlIc, T)i? is to load a band insulator with 
three \g, i) atoms per site on the long lattice and then slowly 
ramp up the short lattice with a bias so that it is energetically 
favorable to have two atoms in the left well and one in the 
right well. Next one can change the state of the right atom 
from |g, to |e, 1)^ by applying a tt pulse of polar- 
ized light resonant with this single-atom transition. The left 
well will be unaffected because the spectrum is modified by 
the interactions (if interactions alone do not provide the de- 
sired selectivity, one could, for example, change the bias of 
the e-lattice). The next step is to change the state of the 
left well from two \g, 1)l atoms populating the lowest two 
vibrational states to \e, I) L\g, I) l both populating the low- 
est vibrational state. This can be accomplished by using tt- 
polarized traveling wave laser light to apply a tt pulse reso- 
nant with the transition between these two many -body states 
il06ll . This results in |e, DlIs, i)L|e, T)i?- One can then 
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FIG. S5: Proof-of-principle experiment to probe RKKY 
interactions in an array of isolated double wells, a, 

Schematic representation of the initial state 
75(1.9, Dl + 1.9, i)fl)|e, i)L|e, T)i?. b, In the limit 
\Vex I <C Jg, the Neel order parameter for the e atoms 

7Ve.(t) « -i COS (%1 ) - I COS (ii^ - Ig^) , Which is 

shown in red for Vex = It exhibits fast oscillations 

with frequency ~ Vex, modulated by an envelope of 
frequency ~ Vj^ / Jg induced by RKKY interactions 

(-3 - 3 COS f 1 shown in blue). 



temporarily shift the g and e lattices relative to each other 
to set U^g interactions to zero, then make Jg nonzero, and 
wait until the g atom evolves into the desired superposition 
73(15, 1)l + 1,9, 1)r) via tunneling. This yields the desired 
state -^ilg, Dl + |.9, i>fl)|e, i)L|e, ^B- 
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